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Abstract 

In 1983 to 84, we proposed a new method in which all stochastic processes concerned are exactly taken into 
account as for range fluctuation of high energy muons, without introducing any approximation which may 
distort fluctuation effect. Now, we call it tentatively Time Sequential Procedure. In 1991, Lipari and Stanev 
proposed another method for the same problem. They divided the problem into two part, soft part and 
hard part, and applied Monte Carlo technique to the latter part. Here, we call their method tentatively Vcut 
Procedure. It is well known that the procedure has been widely utilized in the analysis of high energy muon 
<~i ■ events in KM3 detectors. In the present paper, we examine the limit for applicability to Vcut Procedure 
for the range fluctuation, comparing with Time Sequential Procedure. It is concluded that Vcut Procedure 
^ I give not so different values from Time Sequential Procedure on the survival probability for high energy 
, muons, but it connotes critical problems related to energy measurements on high energy muon events due 
1-^^ • to Cherenkov light signals, owing to the inconsistency of the procedure involved. Thus, we try to the revival 
of Time Sequential Procedure for the establishment of more reliable theory for the observation of Cherenkov 
light whose origins are high and extremely high energy muons. 

Keywords: High energy muon, Muon range fluctuation, Muon energy loss, Muon propagation simulation 



1. Introduction 

The studies on depth intensity problem of high 
energy muons have been one of the most important 
subjects in traditional cosmic ray physics, relating 
to elucidation on unknown character of high energy 
cosmic ray muons, and still have never been lost its 
importance. Menon and Ramana Murthy [1] wrote 
an excellent review on this subject at greater depths 
which describes from the first experiment carried 
out Simizu tunnel in Japan (1940-1945) [2] to the 
last experiment related to the neutrino carried out 
in Kolar Gold Mine, India (Up to 1964) [3]. Bugaev 
et al. [4] have discussed this problem related to 
the charm production mechanism, adding new data, 
from, such as, DUMAND, Baikal, MACRO, LVD, 
NESTOR and others up to 1977. Now, the gigantic 
projects for neutrino astrophysics called as KM3 
detectors are now being developed in the lake [5], 
Antarctic [6] and Ocean [7,8]. In the analysis of 
these KM3 detectors, the depth intensity relation 



for high energy muon is utilized for the confirmation 
of their experimental reliabilities related to other 
experiments different depths [5-8]. 

The theories of range fluctuation of high energy 
muons are indispensable means in the analysis of 
the depth intensity relation of high energy muons 
at certain depths. The theories of range fluctua- 
tion arc studied in three different manners. The 
flrst one is analytical manner [9-19], the second one 
is numerical one [20] and the third one the Monte 
Carlo manner [21-34]. 

At the same time, the theories of range fluctua- 
tion of high energy muons offer essential tools for 
energy determination of high energy muon events 
from neutrino interactions in KM3 detectors. Their 
energies, in particular higher energies, are esti- 
mated from Cherenkov light signals from muon in- 
duced electromagnetic cascade showers, not rather 
than from muons themselves, but, in spite of this 
situation, examination of behaviors of muons them- 
selves are essentially important, because they are 
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origins of electromagnetic cascade showers from dif- 
ferent modes of interactions due to muons. 

At present time, studies on the fluctuation of high 
energy muons have been made by the Monte Carlo 
method, using electronic computers with great per- 
formance, because only this method can clarify fluc- 
tuation characters of muons correctly, while ana- 
lytical method and numerical one provide essen- 
tially their average behavior. The detailed studies 
around fluctuation of high energy muon events are 
inevitable, owing to small number of physical events 
concerned in addition to sharp steepness of the par- 
ent neutrino energy spectrum which arc the origin 
of fluctuation. The studies on the range fluctua- 
tion of high energy muons by Monte Carlo method 
had been made even before appearance of electronic 
computer with great performance [21-25]. How- 
ever, then, they were forced to put more simplified, 
even more artificially assumptions on their stochas- 
tic processes concerned for saving both man powers 
and computer ones at the period for computation. 

In 1983 to 84, Takahashi ct al. [26,27] had de- 
veloped a new Monte Carlo technique where every 
stochastic process for high energy muon concerned 
is treated exactly from the stochastic point of view. 
Namely, the interaction points and the energy divi- 
sion due to interactions concerned (bremsstrahlung, 
direct electron pair production and photonuclear in- 
teraction) are exactly treated in stochastic manner. 
In the present paper, we call it tentatively Time 
Sequential Procedure. 

In 1991, Lipari and Stancv [28] developed another 
technique, from a point of the philosophy of Monte 
Carlo method. They put the diffusion equation on 
the fluctuation in the form of differential-integral 
equation and solve it by Monte Carlo technique. 
They divided the part which is the origin of fluc- 
tuation into two parts, namely, hard part (radia- 
tion loss part) and soft part (continuous energy loss 
part). In radiation loss part only, they treat fluc- 
tuation in Monte Carlo way, but in soft part they 
approximate fluctuation as continuous energy loss, 
supposing small effect of fluctuation. This tech- 
nique has been adopted by subsequent authors [29- 
34]. In the present paper, we call it tentatively Vcut 
Procedure. Now, Vcut Procedure [28-34] has been 
extensively utilized in the analysis of muon neutrino 
events KM3 detectors [for example, 6]. 

However, in our opinion, taking into account of 
the fact that Cherenkov light signal due to high en- 
ergy muons mostly come from the muon induced 
electromagnetic cascade showers whose origins ei- 



ther bremsstrahlung or direct electron pair produc- 
tion or photonuclear interaction, not from muons 
themselves in KM3 detector, the energy determina- 
tions of high energy muon events inevitably include 
more ambiguity in the case of Vcut Procedure, com- 
pared with the case of Time Sequential Procedure 
(see, the section 3 and Conclusion and Outlook ). 

In the present paper, wc try to revitalize Time 
Sequential Procedure in 1993 for the more accurate 
analysis of high energy muon events in KM3 de- 
tector, comparing to the results obtained by Vcut 
Procedure which has been well distributed. Here, 
we restrict our concern to the fundamental and its 
application will be reported in subsequent papers. 

2. Fundamental Structure of Time Sequen- 
tial Procedure 

Here, in order to clarify characteristics of Time 
Sequential Procedure in contrast to Vcut Procedure, 
we reproduce our procedure which had already been 
published in 1984 (in Japanese) [27]. 

2.1. The mean free paths for stochastic processes 
and their resultant mean free path 

Behaviors of high energy muons arc stochas- 
tically determined from elementary processes of 
bremsstrahlung [35] , direct electron pair production 
[36] and photonuclear interaction processes [37]. 
Our idea is very simple and fundamental. 

We treat these processes as stochastic ones as ex- 
actly as possible, without introducing any approxi- 
mation in the energy region in which we are inter- 
estecQ. In our program, these stochastic processes 
are prepared as subroutines and, therefore, they are 
easily replaced by the most advanced ones, if nec- 
essary keeping exactness of our logical structure. 

Let us denote, differential cross sections for 
bremsstrahlung, direct electron pair production, 
and photonuclear interaction, ab {E , Ef,) dEf, , 
ad{E,Ed) dEd and an{E,En)dEn , respectively. 
Here, E denotes the energy of muon concerned, £'{,, 
the energy of 7-ray due to bremsstrahlung, Ed, the 
energy of electron pair due to direct electron pair 
production, £"„ , the energy of emitted particles due 
to photonuclear interaction, respectively. Then, the 



^We adopt the minimum energy of muon for observation 
as 1 GeV throughout paper. Namely, the behaviors of muons 
with larger than 1 GeV arc exactly described, following the 
stochastic probability functions concerned. 
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mean free paths for different stochastic processes 
are energy dependent of the muon concerned and 
they are given as fohows: For brcmsstrahlung pro- 
cesses, 



h (E) = 



1 



A JE^ 



(76 [E, Eb) dEb 
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Here, Eb.min is taken in such a way that Eb,min/ E 
is sufficiently smaller than Emin/E. Here, Emin de- 
notes the minimum energy of the muon for obser- 
vation. Emin is taken as f GeV throughout present 
papei0. The integrations for direct electron pair 
production and photonuclear interaction are per- 
formed over kinematically allowable ranges. For di- 
rect electron pair production processes, 
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For photonuclear interaction processes, 
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, where A'^ and A denote Avogadro number and 
atomic mass number, respectively. 

Also, Atotai {E) , the resultant mean free path for 
these stochastic processes are given as. 
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2.2. Determination of the kind of the stochastic 
processes and the real free path for the stochas- 
tic processes concerned 

By using Eq. (f) to (4), we can determine the 
interaction points of muons for different stochastic 
processes in the following. The first, for the pur- 
pose, lets us define {E) and {E) as follows; 
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^Such the limitation is absolutely necessary for the exact 
treatment 7-rays due to brcmsstrahlung in the energy region 
where we are interested owing to the infrared catastrophic 
characteristics of brcmsstrahlung. In the present paper, the 
behavior of muons with energies larger than 1 GeV related 
to brcmsstrahlung are exactly treated in stochastic manner, 
as well as those for direct electron pair production and pho- 
tonuclear interaction. 



The second, we sample randomly ^1 , a uniform 
random number between (0,1). If ^1 < (E) , 
then we recognize the interaction occurs due to 
brcmsstrahlung. If S,b (E) < Ci < (E) , then 
we understand the direct electron pair production 
occurs. If ^1 > ^d{E), then, we understand that 
photonuclear interaction occurs. 

Again, we sample a new ^2 , randomly from uni- 
form random number between (0,1). Then, we 
can determine the interaction points At for the 
specified stochastic processes according the follow- 
ing criterion. In the case of the occurrence of 
brcmsstrahlung processes ( for ^1 < (E) ), 



Atb = -Xtotai (E) log£,2 



(7) 



In the case of the occurrence of direct electron pair 
production processes (for ^j, [E) < ^1 < {E) ), 



At. 



'Xtotal (E) log^2 



(8) 



In the case of the occurrence of photonuclear inter- 
action processes ( for ^1 > (E) ) 



At„ = -Xtotai (E) log^2 



(9) 



We are taken into account of ionization loss addi- 
tion to stochastic energy losses, by using Tamura's 
method [38]. 

2.3. Determination of the emitted energy loss (Eb 
or Ed or En) from the specified stochastic pro- 
cesses 

Under the determination of the interaction points 
due to the specified stochastic processes, by using ^1 
and ^2, in the previous subsection, here, the energy 
losses Eb , Ed and En from the specified stochas- 
tic processes are given as follows. For sampled ^3 
which is obtained randomly from the uniform ran- 
dom number between (0,1), we solve the following 
equations for respective interactions in order to ob- 
tain Eb or Ed or En- For brcmsstrahlung process. 
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For photonuclear interaction, 
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(12) 
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In Eq. (10) to (12) for respective interaction, the 
quantities to be obtained are Eb^d,n, emitted en- 
ergy losses for given i?, energy of the muon con- 
cerned. The quantities of Eh_d,n are extended to 
Eb,d,n,min to Eb,d,n,max- For Sampled ^3 we can 
solve these equation numerically and obtain Ei,^d,n 
finally. Thus, we can determine Eb_d,m the energy 
losses for the specified stochastic process at deter- 
mined interaction point. A flow chart for the fun- 
damental structure of Time Sequential Procedure is 
given Fig.l. 
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Figure 1; Flow Chart for the fundamental structure of Time 
Sequential Procedure. 



2.4- On the validity of Time Sequential Procedure 

Generally speaking, the verification of the valid- 
ity of the Monte Carlo method concerned is pretty 
difficult. For the purpose, it is desirable to compare 
the physical results obtained by Time Sequential 
Procedure with the corresponding results obtained 



by the analytical method which is methodologically 
independent of the Monte Carlo method concerned. 

Once, Misaki and Nishimura [17] had developed 
an analytical theory for range fiuctuation of high 
energy cosmic ray muons based on the Nishimura- 
Kamata formalism on electron shower theory to ap- 
ply to study depth intensity relation muon under- 
ground. Takahashi et al. [26,27] had compared the 
depth dependence of the average energies of muons 
under the incident energy spectrum of the muons on 
the surface with the different indices (7= 2, 3, and 
4 ) obtained the analytical theory [17] with the cor- 
responding quantities obtained by Time Sequential 
Procedure. 




Depth(g/cm 



Figure 2: The average energies of the muons. The lines 
denote Misaki and Nishimura, while symbols ours. 



We reproduce them in Fig. 2 from the previous 
work [27]. The agreement between them is quite 
well, taking into account of the difference in both 
the cross sections concerned and their numerical 
evaluation method. Namely, we can say the logi- 
cal structure of Time Sequential ProcedureOn the 
validity of is well established from the point of the 
validity of Monte Carlo method concerned. 

2.5. Examples directly derived by the Time Sequen- 
tial Procedure 
2.5.1. The diversity of individual muon behavior 
'Needle' structure of the energy losses from 
high energy muons 
In subsections 2.2 and 2.3, by random sampling 
procedure, we show how to determine both inter- 
action points for the specified stochastic processes 
and subsequent their energy losses for the inter- 
actions concerned. In the present subsections, we 
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show some examples of the energy losses along the 
passage of the high energy muons. 

In Fig. 3 to Fig. 13, we show three different cate- 
gories of the 'Needle' structure of the energy losses 
for the same primary energies. 'Needles' denote 
fractional energy losses due to a specified stochas- 
tic processed such as, bremsstrahlung, direct elec- 
tron pair production or photonuclear interaction at 
respective interaction points. We illustrate several 
typical structures of the energy loss of muons which 
start, having primary energy Eq and reaching Emm, 
namely, the behaviors with the shortest range, with 
the longest range and with the average-like range 
for the same primary energy and the same starting 
point. The shortest range denotes the muon with 
the shortest range among all sampled muons, while 
the longest range does the muon with the longest 
range among all sampled muons, and the average- 
like range does the muon with the range whose is 
the nearest to the averaged range among all sam- 
pled muons. Total sampling numbers per respective 
primary energy are 100,000. 

In these figures, we can recognize the diversi- 
ties of muon behaviors for the same primary muon 
energies with regard to their ranges (or their en- 
ergy losses). All interaction points due to the pro- 
cesses of bremsstrahlung, direct electron pair pro- 
duction and photonuclear interaction and all energy 
losses due to these elementary processes at respec- 
tive points due to these processes are recorded. In 
order to clarify the diversities among the real range 
distributions (or real energy loss distributions), we 
examine the muons with the shortest range, the 
muons with the longest range and the muons with 
the average-like range in more detail. 

In Table 1, we show numerically the characteris- 
tics of an individual muon with the shortest range, 
the average-like range, the longest range, in addi- 
tion to their average range for IG^^eV, lO^^eV and 
lO^^eV. 

In Fig. 3 to Fig. 5, we give the characteristic be- 
haviors with the shortest range, the average-like 
range and the longest range, that is, their energy 
loss for the specified interaction as the function 
of the depth traversed for the primary energy of 
lO^^eV in watefl In these figures, we utilize the 
same scale in depth to clarify the diverse behav- 



■^In order to understand the situation visually the char- 
acteristic behaviors of high energy muons which are shown 
Fig. 3 to Fig. 13, we suggest the readers to look at the pictures 
with colors in the WEB page. 



iors by the same incident energies, namely, those 
with the shortest range, with the average-like range 
and with the longest range, respectively. In figures, 
the abscissa denotes the depths where the specified 
interactions occur. The 'needles' (expressed in ordi- 
nate) with different colors at different depths denote 
ratios of the energy losses due to direct electron pair 
production (green, d), bremsstrahlung (red, b) and 
photonuclear interaction (blue, n) to their primary 
energy, respectively. The abrupt changes in them 
are due to the catastrophic energy losses for muons 
(see, footnote 2). It is easily understood that one 
sees the fluctuation effect rather week in the energy 
of IQi^eV. 

It is seen from figures and Table 1 that there 
is not so big difference between the case with 
the shortest range and one with the longest range 
for lO^^eV. In the case with the shortest range 
(Fig. 3), we find two catastrophic energy losses (at 
~ 910 meters and ~ 1870 meters) due to two 
bremsstrahlungs play the important role in the 
range. In the case with the average-like range 
(Fig. 4), we can find one catastrophic energy loss 
due to bremsstrahlung at ^ 1.48 kilometer. How- 
ever, in the case of the longest range (Fig. 5), we 
cannot find the catastrophic energy losses due to 
bremsstrahlung and, instead, we can find that al- 
most energy losses are due to many number(^ 300) 
of direct electron pair production events. 

In Fig. 7 to Fig. 9, we give the typical diversi- 
ties for primary energy of 10 ^^eV similarly for pri- 
mary energy of lO^^eV. In these figures, the diver- 
sities for the shortest range, the average-like range 
and the longest range are compared explicitly ex- 
pressed in the same scale. Fig. 6 shows the same 
in Fig. 7 in extended scale. Combined with Table 
1, the shortest range, ^ 940 meter (Fig. 6), is far 
shorter compared with the longest range, ~ 35.0 
kilometers (Fig. 9). It is seen from Fig. 6 and the 
Table 1 that bremsstrahlung plays a decisive role 
as the cause of catastrophic energy loss in the case 
of the shortest range, ( ^ 96.5% of total energy up 
to^ 450 meters). 86.6% of the total energy is lost 
by 2 bremsstrahlungs, 13.4% by 367 direct electron 
pair productions and 5.9 x 10""^% by 1 photonu- 
clear interaction. In Fig. 9, we give the case for 
the longest range. Here, large numbers of direct 
electron pair production with rather small energy 
loss play an important role, as shown similarly in 
Fig. 5. Here, 80.2% of the total energy is lost by 
13722 direct electron pair productions, 7.53% by 71 
bremsstrahlungs and 11.1% by 105 photonuclear in- 
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Figure 3: The fractional energy loss with the shortest range as the function of the depth for lO^^eV muon. A line graph in the 
upper denotes fractional muon energy as the function of the depth, [b] denotes the fractional energy loss due to bremsstrahlung, 
[d] due to direct electron pair production and [n] due to photonuclcar interaction, [b], [d] and [n] are utilized as the same 
meaning up to Figure 13. 
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Figure 4: The fractional energy loss with the average-like range for lO^^eV muon muon together with the fractional energy of 
the muon. 





1.5 2 2.5 

Deptn in water(km) 



Figure 5; The fractional energy loss with the longest range for 10 eV muon muon together with the fractional energy of the 
muon. 



teractions. In Fig. 8, combined with Table 1, we give 
the case with the average-like range. Here, 23.5% 
of the total energy is lost by 5489 direct electron 
pair productions, 75.0% by 49 bremsstrahlungs and 
0.93% by 37 photonuclear interactions, while in the 
real averages (100,000 samples), 47.4% of the total 
energy is lost by 6800 direct electron pair produc- 
tions, 35.3% by 48.1 bremsstrahlungs and 16.7% by 
55.0 photonuclear interactions. 

In Fig. 10 to Fig. 13, combined with Table 1, we 
show the similar relations for lO^^eV muons as 
shown in lO^^eV. We can say the case with the 
shortest range in Fig. 10 (or Fig. 11) has a strong 



contrast to that with the longest range. The manner 
of the energy loss in Fig. 10 is drastic with two big 
catastrophic energy losses due to bremsstrahlungs 
(^ 0.8km and 4.74km), while that in Fig. 13 is very 
moderate with no catastrophic energy loss. The 
shortest range, ^ 7.7 kilometers (Fig. 11), is far 
shorter compared with the longest range, ^ 57.8 
kilometers (Fig. 13). It is seen from Fig. 10 and 
Table 1 in the case of the shortest range that 
bremsstrahlung plays a decisive role as the cause 
of catastrophic energy loss. 87.5% of the total en- 
ergy is lost by 28 bremsstrahlungs, 11.9% by 5760 
direct electron pair productions and 0.623% by 40 
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Table 1: The details of the characteristics on the muons with the shortest range, the average-like range, the longest range and 
the average range. 



Eo = lO^^eV 


Range 
[km] 


Energy loss 
by brems[eV] 


Number of 
interaction 


Energy loss by 
direct pair[eV] 


Number of 
interaction 


Energy loss 
by nuclear [eV] 


Number of 
interaction 


<Average> 


2.43 


1.10x10^1 


4.74 


1.57x10^1 


243 


4.54x10"' 


3.44 


Average-like 


2.43 


8.97x10"' 


4 


1.34x10" 


221 


8.86x10"' 


3 


Shortest 


2.09 


2.15x10" 


6 


1.52x10" 


208 


3.72x10"' 


3 


Longest 


3.14 


3.80x10" 


5 


1.04x10" 


299 


1.19x10'^ 


3 


Eo = lO^^eV 




<Average> 


1.78x10^ 


3.53 xlO^"* 


48.1 


4.74x10^^ 


6.80x10-^ 


1.67x10" 


55 


Average-like 


1.78x10^ 


7.50x10^'' 


49 


2.35x10" 


5489 


9.31x10^^ 


37 


Shortest 


9.44x10"^ 


8.66x10^" 


2 


1.34x10^^ 


367 


5.90x10"' 


1 


Longest 


3.50x101 


7.53x10"^ 


71 


8.02x10^^ 


13722 


1.11x10" 


105 


Eo = lO^^eV 




<Average> 


3.28x10^ 


3.37x10"'' 


108 


4.39x10"'' 


2.57x10* 


2.25x10"' 


172 


Average-like 


3.28x10^ 


1.68x10"'' 


118 


5.58x10"'' 


29321 


2.74x10"' 


196 


Shortest 


7.72x10" 


8.75x10"'' 


28 


1.19x10"' 


5760 


6.23x10"' 


40 


Longest 


5.78x10^ 


5.71x10^" 


162 


5.77x10" 


46542 


3.66x10" 


277 



photonuclear interactions. In Fig. 13, we give the 
case with the longest range. Here, 57.7% of the 
total energy is lost by 46542 direct electron pair 
productions, 36.6% by 277 photonuclear interac- 
tions and only 5.71% by 162 bremsstrahlungs in 
the complete absence of catastrophic energy losses. 
In Fig. 12, we give the case with the average-like 
range. Here, 55.8% of the total energy is lost 
by 29321 direct electron pair productions, 16.8% 
by 118 bremsstrahlungs and 27.4% by 196 pho- 
tonuclear interactions, while, in the real averages 
(100,000 samples), 43.9% of the total energy is lost 
by 2.57x lO'^ direct electron pair productions, 33.7% 
by 108 bremsstrahlungs and 22.5% by 172 photonu- 
clear interactions. Thus, it can be concluded that 
the diversity among muon propagation with the 
same primary energy should be noticed. 

2.5.2. Average characteristics of high energy muons 
with the shortest range, the average-like 
range, the longest range around the average 
range 

In Table 2 (a), we give the ratios of energy 
losses due to respective stochastic processes to to- 
tal energy loss in the typical ranges {the average, 
the average-like, the shortest and the longest) for 
lO^^eV, lO^^eV and lO^^eV. It is clear from the ta- 
ble that, averagely speaking, high energy muons are 
lost ~ 50% in the direct electron pair production, 
^ 30% in bremsstrahlung and ~ 20% in photonu- 



clear interaction. It is clear from Table 2(a) that the 
muon with the longest range loses ~ 70 % in direct 
electron pair productions, in long chain of electron 
pairs induced electromagnetic cascade showers with 
rather smaller energies, while the muon with the 
shortest range loses ~ 70% to 100% of their energy 
in a few number of bremsstrahlungs (catastrophic 
energy loss). 

In Table 2(b), the ratios of fractional energy lose 
for specified stochastic processes are divided by 
the corresponding averaged ones. It is clear from 
the Tabic 2(b) that the divisions of energy loss to 
the specified processes in the average-like range arc 
clearly different from that of the averaged. It shows 
that the energy divisions for different processes arc 
different, even if the ranges are same. This fact 
makes the ejection of the Cherenkov light influence, 
even if the paths of the high energy muons are same. 

In Table 2(b), it is also clear from the character- 
istics of the typical showers from the point of en- 
ergy dissipation that energy losses ratios of showers 
concerned to their averages due to bremsstrahlung 
in the shortest ranges lose their energies arc 2.96, 
2.53, 2.17 for primaries lO^^cV, lO^^eV and lO^^cV, 
respectively. Namely, these showers with the short- 
est range essentially lose their energies almost due 
to bremsstrahlung, while the corresponding ratios 
due to direct electron pair production in the show- 
ers with the longest range are 1.33, 1.52, 1.49. 
Also, these showers with the longest range lose their 
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Figure 6: The fractional energy loss with the shortest range for lO^^cV muon together with the fractional energy of the muon. 
The figure is a magnification of Figure 7. 
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Figure 7: The fractional energy loss with the shortest range for lO-'^^eV muon together with the fractional energy of the muon. 
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Figure 8: The fractional energy loss with the average-like range for lO^^cV muon together with the muon energy. 
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Figure 9: The fractional energy loss with the longest range for lO^^eV muon together with the fractional energy of the muon. 



pretty energies owing to direct electron pair produc- 
tion. 



2.5.3. Range Distributions and Hypothetical Range 
Distributions for high energy muons 
As shown, for example, in Fig. 3 to Fig. 13, we 
can pursue three kinds of the typical types of the 
behaviors of high energy muons with definite pri- 
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Figure 10: The fractional energy loss with the shortest range for 10^*cV together with fractional energy of the muon. The 
figure is a magnification of Figure 11. 
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Figure 11: The fractional energy loss with the shortest range for 10^** eV muon together with the fractional energy of the muon. 
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Figure 12: The fractional energy loss with the average-like range for 10 eV muon together with the fractional energy of the 
muon. 
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Figure 13: The fractional energy loss with the longest range for 10^'^eV muon together with the fractional energy of the muon. 




mary energies in stochastic manner exactly, record- 
ing the locations of the interaction points for spec- 
ified interactions and their dissipated energies ex- 
actly. However, we pursue the behaviors of all sam- 



pled muons exactly, including three different types 
of the muons. We can construct the range distri- 
butions from ensemble of 100,000 individual muons 
for respective primary muon's energy, as shown in 
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Table 2: The ratios of energies transferred from bremsstrahlung, direct electron pair production and photonuclear interaction 
to the total energy loss (a) and their ratios expressed in respective average values (b). 





Brems 


Direct Pair 


Nuclear 


J-JQ — -1- U C V 


(a) 


(b) 


(a) 


(b) 


(a) 


(b) 




3.37x10-1 


1.00 


5.26x10-1 


1.00 


1.37x10-1 


1.00 


A VPrFI ffp-lilcp 


2.87x10-1 


0.872 


4.30x10-1 


0.688 


2.83x10-1 


2.51 


Shortest 


5.32x10-1 


2.96 


3.76x10-1 


9.14 X 10-1 


9.20x10-^ 


2.96 X 10-^ 


Longest 


3.50x10-^ 


0.875 


9.54x10-1 


1.33 


1.10x10-^ 


5.87 X 10-^ 


E„ = IQ^-'eV 




<Average> 


3.40x10-1 


1.00 


4.98x10-1 


1.00 


1.62x10-1 


1.00 


Average- like 


7.54x10-1 


1.25 


2.37x10-1 


0.960 


9.36x10-'' 


0.600 


Shortest 


8.66x10-1 


2.53 


1.34x10-1 


0.178 


5.90x10-^ 


0.321 


Longest 


7.62x10-^ 


0..541 


8.11x10-1 


1.52 


1.12x10-1 


0.364 


Eo = lO^^eV 




<Averagc> 


3.24x10-1 


1.00 


4.59x10-1 


1.00 


2.17x10-1 


1.00 


Average- Uke 


1.68x10-1 


1.30 


5.58x10-1 


0.858 


2.74x10-1 


0.848 


Shortest 


8.75x10-1 


2.17 


1.19x10-1 


0.634 


6.23x10-'' 


2.20 X 10-^ 


Longest 


5.71x10-^ 


0.209 


5.77x10-1 


1.49 


3.66x10-1 


1.13 



Fig. 14. In the figure, we give P{R;Eq), the proba- 
bilities for the range distribution in water with pri- 
mary energies, lO^^eV to lO^^eV and lO^^eV in wa- 
ter whose minimum energy is 10®eV(lGeV), respec- 
tively. It is clear from the figure that the width of 
the range distribution increases rapidly, as their pri- 
mary energy increases. Also, as the primary energy 
decreases, the width of range distribution becomes 
narrower and approaches to a 5 function-type, the 
limit of which denotes no fluctuation. It is interest- 
ing that the range distributions can be well approx- 
imated as the normal distribution above ~ lO^^eV 
where the total Cherenkov light yields comes al- 
most from the muon induced electromagnetic cas- 
cade showers and they are given as, 

, where Eq, R, < R > and a are primary energy, 
real range, the average value of ranges and the stan- 
dard deviations, respectively. Their average ranges, 
standard deviations and relative variances (stan- 
dard deviations divided by averages) in water are 
given in Table 3. Also, it is interesting that their 
relative variances decrease slightly as their primary 
energies increase. It should be noticed from Table 
3 that the standard deviation increases as primary 
energy increases, but, the relative variance of the 
range distribution decreases inversely. 

In order to examine each characteristic of the 
stochastic process, such as the bremsstrahlung, di- 
rect electron pair production and photonuclear in- 




Depth(km) 

Figure 14: Range distributions for IQi'^eV to IQl^eV muons 
in water. The minimum observation energies are taken as 
10^ eV. Each sampling number is 100,000. 



Table 3: The average values, the standard deviations and the 
relative variances of the range distributions of muons from 
IQlleV to IQl^eV in water. 



Eo [cV] 


<R> [km] 


a [km] 


a/ <R> 


IQii 


3.56x10-1 


2.52x10-^ 


7.07x10-^ 


IQi^ 


2.43 


4.71x10-1 


1.94x10-1 


IQl'' 


7.28 


2.02 


2.78x10-1 


IQi'i 


1.26x101 


3.49 


2.77x10-1 


10^^ 


1.78x101 


4.57 


2.57x10-1 


10l« 


2.30x101 


5.41 


2.36x10-1 


IQi'^ 


2.79x101 


6.14 


2.20x10-1 


10i« 


3.29x101 


6.81 


2.07x10-1 



10 



teraction, wc construct the hypothetical range dis- 
tribution in which a specified stochastic process 
only is assumed to occur and the other two stochas- 
tic processes arc assumed not to occur. To clarify 
the characteristics of the specified stochastic pro- 
cesses, we can compare this hypothetical range dis- 
tribution with that of real range distribution in 
which every specified stochastic process is realized 
as the competition effect among these three pro- 
cesses. We compare the hypothetical range distri- 
bution with the real range distribution in Fig. 15 to 
Fig.l7. 

In Fig. 15, we compare three different hypothet- 
ical range distributions with the real range dis- 
tribution for primary energy of 10^'^eV. Here, the 
symbol d in these figures means the hypotheti- 
cal range distribution in which only direct electron 
pair production is taken into account and both the 
bremsstrahlung and photo nuclear interaction are 
neglected. The symbols b and n have similar mean- 
ing to that of d. The symbol t means the real range 
distribution in which all interactions are taken into 
account (The true distribution). From the shapes 
of the distributions and their maximum frequencies 
for different stochastic processes in the figures, it is 
clear that energy losses in the direct electron pair 
production arc of small fiuctuation, while both the 
bremsstrahlung and photonuclear interaction are of 
bigger fiuctuation, and the fluctuation in photonu- 
clear interaction becomes bigger when compared 
with bremsstrahlung as primary energy increases. 
The smaller fluctuation in direct electron pair pro- 
duction suggests us that energy loss from this pro- 
cess may be treated as something like continuous 
energy loss in the special situation 

2.5.4- Other physical quantities obtained from 
Time Sequential Procedure 

In Fig. 18 to Fig. 20, we give the survival proba- 
bilities for different cutoff energies with primary en- 
ergies of lO^^eV, lO^^eV and lO^^eV, respectively. 
The values for cutoff energies are given in respective 
figures. The sampling number utilized is 100,000 for 
each primary energy. 

In Fig. 21 to Fig. 23, we give the differential 
energy spectrum of muons for primary energies, 
lO^^eV, lO^^eV and lO^^eV, respectively. The en- 
ergy spectrum of the survival muon obtained by 
Vcut Procedure is surmised to be different from 
that by Time Sequential Procedure, in particular, 
in higher energies. Their magnitude is smaller than 
that by Time Sequential Procedure and their shapes 
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Figure 15: Hypothetical range distributions in water for 
lO^^eV muons together with the real range distribution. 
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Figure 16: Hypothetical range distributions in water for 
10^^ eV muons together with the real range distribution. 
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Figure 17: Hypothetical range distributions in water for 
lO^^eV muons together with the real range distribution. 



are different from that by us. The reason for the 
difference is discussed briefly in the following sec- 
tion. 
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Figure 18: The survival probabilities for lO^'^eV muon. 
Curves labels correspod to following set of cutoff energies: 
(l)109cV, {2)10l0cV, (3)10"cV. 
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Figure 19: The survival probabilities for lO^^cV muon. 
Curves labels correspod to following sot of cutoff energies: 
from (l)10''eV to (6)10"eV. 
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Figure 20: The survival probabilities for lO^^eV muon. 
Curves labels correspod to following set of cutoff energies: 
from (l)109eV to {9)10"eV. 



Figure 21: Energy spectrum in water at the different depths, 
initiated by 10^'^cV muons. 
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Figure 22: Energy spectrum in water at the different depths, 
initiated by lO^^cV muons. 
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Figure 23: Energy spectrum in water at the different depths, 
initiated by lO^^eV muons. 



3. Vcut Procedure: The fundamental struc- 
ture and points at issue 

3.1. Vcut Procedure as an approximated Time Se- 
quential Procedure 
Nobody has verified tlie validity of the logical 
12 structure of Vcut Procedure from the view point 



of the fundamental idea of Monte Carlo method. 
The plausible reason why Vcut Procedure seems to 
give right results is simply that many people have 
utilized this teehnique and they obtain essentially 
same results. However, this does not always guar- 
antee that Vcut Procedure is logieally correet in the 
sense of Monte Carlo method, because it is quite 
natural for different authors to reach same results, 
if they utilize the same principle which, in same 
case, may be wrong. That is, the same results ob- 
tained by the Monte Carlo method due to different 
authors should be carefully examined in the spirit 
of Monte Carlo method. 

As similarly in our examination (see, section 2.4), 
the validity of Vcut Procedure in the sense of Monte 
Carlo method should be examined by methodologi- 
cally independent of Monte Carlo method, namely, 
by either analytical method or numerical one, focus- 
ing on the average values of several physical quan- 
titiefl It is rather easy to develop a numerical 
theory based on the same principle and compare 
the results obtained by the numerical theory with 
corresponding results obtained by Vcut Procedure. 
Here, we examine the fundamental structure of Vcut 
Procedure from a different point of view. 

Lipari and Stanev [28] and subsequent au- 
thors, P.Antonioli, S.Iyer et al., Klimushin et al., 
D.Chirkin et al., S.Bottai et al. [29-34] formulate 
Vcut Procedure as follows: 



dE 
dx 



dE 
dx 

N 



+ 



soft 



dE 



dx 



hard 



a{v,E) 
dv ■ V- 



A 



E 



dv ■ 



dv 
da {v, E) 
dv 



(14) 



, where v denotes the fractional emitted energy. 
They introduce Vcut , a certain constant value, into 
the diffusion equation, in such a way that the effec- 
tive energy loss, for example, 7-ray above Vcut x E, 
is treated in stochastically only hard part, while 



''For example, in order to confirm the validity of the one 
dimensional Monte Carlo method for electromagnetic cas- 
cade shower, we should compare the averaged transition 
curves and/or track length for shower particles under Ap- 
pro. A with the corresponding results obtained by the ana- 
lytical theory. In three dimensional cases, the mean square 
of angular and lateral spreads for shower particles should be 
compared with each other for the validity of the Monte Carlo 
method [39-42]. 



that below Vcut x _B is put into continuous energy 
loss {soft part), which is simply subtracted from the 
muons concerned. 

Here, let us summarize the values of Vcut utilized 
in Vcut Procedure in the following, [a] Lipari and 
Stanev adopt Vcut— 0.01 [28], [b] Antonioli et al. 
adopt Vcut ~ 10^'^ [29], [c] Dutta et al. adopt Vcut 
=10-3 [31], [d] Sokalski et al. adopt Vcut= 10"^ to 
0.2 [32], [e] Chirikin and Rohde adopt Vcut = 10"" 
to 10-3 [32]. 

The problems to be involved into Vcut Procedure 
are as follows: 

[A] It should be pointed that the separation of 
soft part from hard part is treated in inconsistent 
manner in Vcut Procedure as for fixed energy muon. 
Namely, the muon with the some energy is treated 
in soft part in some case, while the muon with the 
same energy is treated in hard part in another case. 
Such the treatment lack in consistency for descrip- 
tion on muon behavior, because the effectiveness of 
fluctuation depends on the absolute values of muon 
energies. 

Vcut Procedure pursues the change of energy state 
by step by step method with regard to the depth dx. 
Consequently, by the constancy of Vcut (10-"* un- 
der examination), their stochastic energy loss part 
[hard part) shifts toward lower energy region, as 
dx advances. In other words, as already mention, 
the muons with some energy belongs to hard part 
(stochastic energy loss part) at certain depth, but 
belongs to soft part (continuous energy loss part) at 
another depth, owing to the shift of the boundary 
line between hard part and soft part. Such a de- 
scription on the behavior of the muon in Vcut Pro- 
cedure clearly lack in consistency. 

For example, comparing Fig. 24(a) with 
Fig. 24(b), it is clear that the region from 10^'^eV 
to 10"eV for Eq = W^eV belongs to soft part, 
while the same region belong to hard part for 
Eo = lO^^eV. This is also an example that the 
stochastic process is not treated in the unified 
manner. 

On the contrast to Vcut Procedure, there are 
no separation of soft part from hard part in 
Time Sequential Procedure and only hard part ex- 
ists. Namely, all plausible processes are treated in 
stochastic manner exactly in this procedure. 

[B] For higher energy muons and/or larger Vcuti 
for example, in the case of bremsstrahlung, the 
emitted higher 7-rays may be contained in soft part 
in which they are subtracted from the muon con- 
cerned, being treated continuous energy loss and 
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Figure 24: (a) The separation of stochastic energy loss part 
(hard part) from Continuous energy loss part (soft part) in 
the case Vcut = 10~* for 10^* eV. (b) The separation of 
stochastic energy loss part (hard part) from Continuous en- 
ergy loss part (soft part) in the case Vcut = 10"'* for lO^^eV 



don't contribute to the muon's future behavior any 
more. However, some part of such the higher 7- 
rays may not be consumed as dissipated energies, 
if Vcut sets up more smaller. Then the muons con- 
cerned should maintain them in hard part and their 
energy loss may be treated in stochastic manner. 
As the result of it, the emitted 7-rays may be cor- 
rectly taken into account in hard part so that the 
muon concerned can maintain higher energy than 
that in the case of larger Vcut and, consequently, 
more muons may survive than in the case of larger 

Vcut- 

Namely, the present situation with larger Vcut in 
higher energies may induce smaller survival proba- 
bility for the muon concerned. Thus, in the case of 
higher energy muons and/or larger Vcut, we expect 
the deformed high energy 7-rays (or electrons) spec- 
trum, which, in turn, may result in the deformed 
Cherenkov light spectrum, compared with the most 
probable spectrum which can be obtained by Time 
sequential Procedure. This is inevitable problem, 
because these high energy 7-rays (electrons) spec- 
trum are exclusively the sources of Cherenkov light 
in higher energies. (See, Conclusion and Outlook). 

3.2. Comparison of the survival probabilities ob- 
tained by Time Sequential Procedure with those 
obtained by Vcut Procedure 

In the previous section, we examine the incon- 
sistency problems involved in Vcut Procedure. How- 
ever, Vcut Procedure may be useful within their lim- 
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Figure 25: The comparison of our result with that of Lipari 
and Stanev [28ll . The survival probabilities of muons of en- 
ergy from 1 TeV to 10® TeV. The numerical figures attached 
each curve denote the primary energies. Curves labels corre- 
spond to following set of primay energies of muon: (l)lTeV, 
(2)10TeV, (3)102TeV, (4)10^TeV, (5)10-*TeV, (6)105TeV, 
(7) 10® TeV. Symbols are due to Lipari and Stancv and curves 
are due to ours. 
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Figure 26: The comparison of our results with that of 
Klimushin et al[? ]. The continuous lines are obtained by us, 
while symbols are readout from those by Klimushin et al for 
primary energies from 10^^ eV to 3 X 10*® eV. The numeri- 
cal figures attached each curve denote the threshold energy is 
10 TeV. Curves labels correspond to following set of depths: 
(1)1. 15km, (2)3.45km, (3)8. 05km, (4)12.65km, (5)17.25km, 
(6)21. 39km. 



itation, because it may be regarded as an approxi- 
mation of Time Sequential Procedure. 

In Fig. 25 and Fig. 26, we give the comparison of 
our results by Time Sequential Procedure for sur- 
vival probabilities with those of Lipari and Stancv 
and those of Klimushin et al. by Vcut Procedure, 
respectively. We discuss the agreement or disagree- 
ment between the results obtained by Time Sequen- 
tial Procedure and Vcut Procedure. The agreement 
between Lipari and Stanev's (Fig. 25) and ours is 
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quite well in the energies from 1 TeV to lO^'TeV, 
while the disagreement between them become clear 
beyond lO^eV. As indicated in the beginning of [B], 
the 7-rays contained soft part, should be involved 
in the muon spectrum, if the stochastic processes 
concerned are taken into account. Therefore, this 
fact leads increase of survival probability and muon 
spectrum. 

The agreement and disagreement between 
Klimushin et al's (Fig. 26) can be explained simi- 
larly. 

Also, it should be noticed that fluctuation effect 
depend entirely depend on the absolute values of 
muon's energy. 

In Vcut Procedure, Monte Carlo procedure, the 
calculations are performed by step by step method 
in dx, which inevitably introduce uncertainties due 
to the accumulation effect in calculation error. 
One may call their procedure differential method. 
On the contrary to it. Time Sequential Procedure 
can determine interaction point directly (integral 
method). Consequently, its accuracy is indepen- 
dent of the errors coming from the accumulation 
effect coming from dx. One may call our procedure 
as integral method on the contrast to differential 
method. 



4. Conclusion and Outlook 

We have showed the validity of Time Sequential 
Procedure in the Monte Caro method, through the 
comparison of our results with those obtained by 
the analytical theory which is methodologically in- 
dependent. Also, we have examined the limitations 
for application imposed upon Vcut Procedure. 

In KM3 detectors, we catch signals from 
Cherenkov light due to high energy muons. It 
should be emphasized here that Cherenkov light 
signal caught by KM3 comes exclusively muon in- 
duced electromagnetic cascade shower and the con- 
tribution from muon themselves is negligible. 

Now, we discuss this problem in more detail. 
There are three different types of electromagnetic 
cascade showers among Cherenkov light which are 
produced from muon induced electromagnetic cas- 
cade shower. The first one is 7-ray induced electro- 
magnetic cascade showers due to bremsstrahlung, 
the second is e-|- and e- induced electromagnetic 
shower due to direct electron pair production and 
the third is ttq — 27 induced electromagnetic cas- 
cades showers due to pliotonuclear interaction. In 



our procedure, the electromagnetic cascade show- 
ers are also exactly simulated which are origins of 
Cherenkov light. 
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Figure 27: Ratio of Cherenkov lights due to the accompanied 
cascade showers to the Cherenkov hght. 



In Fig. 27, we give the ratios of Cherenkov light 
due to muon induced electromagnetic cascades 
showers to the total Cherenkov light as the func- 
tion of the depth traversed for lO^^eV to 10^^ cV. 
It is clear from the figure that ^ lO^^eV, the most 
of Cherenkov light comes from muon itself, while 
above lO^^eV the most Cherenkov light comes from 
the muon induced electromagnetic cascade showers 
and, consequently, Cherenkov light from the orig- 
inal muon is completely negligible. Just this situ- 
ation around Cherenkov light production strongly 
imposes the limit for application of Vcut Procedure 
to estimation of the energy of muon events in KM3 
detectors. 

Thus, the origins of the electromagnetic cas- 
cade shower arc cither 7-ray (in the case of 
bremsstrahlung and photonuclear interaction) or 
electrons (in the case of direct electron pair pro- 
duction) and, consequently, in order to obtain 
Cherenkov light yield correctly, it is necessary 
to know the behaviors of electromagnetic cascade 
showers due to different stochastic processes which 
are sources for Cherenkov light and, furthermore, 
to know the energy spectrum of 'primary' electron 
and 7-ray which are the origins of electromagnetic 
cascade showers. In Time Sequential Procedure, we 
exactly the muon behaviors in stochastic manner, 
without introducing soft part and. consequently, we 
obtain accurate muon energy spectrum at arbitrary 
depths as well as the energy spectrum of 'primary' 
(in the sense of origin of electromagnetic cascade 
showers) electrons and 7-rays at arbitrary depths. 
However, in Vcut Procedure, one cannot obtain accu- 
rate 'primary' electrons and 7-ray energy spectrum 
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due to strong limitation involved in their technique, 
which is examined in section 3. In particula, the 
emitted energies from soft part are difficult to han- 
dle accurately for higher muon energies and rather 
larger Vcut hi Vcut Procedure. 

Thus, it is rather easily concluded that, compared 
with Time Sequential method, Vcut Procedure in- 
volve the difficulty around accurate muon energy 
spectrum, but also 'primary' 7-rays and electron 
energy spectrum which are sources of Cherenkov 
light, which may distort finally the fluctuation ef- 
fect in Cherenkov light. 

Now, our discussion still remains to single high 
energy muon problem. However, really, muon does 
not exist singly, but they exist in the form of energy 
spectrum which is directly reflection of parent neu- 
trino spectrum. The existence of energy spectrum 
of muons brings more difficulty into the elucidation 
of fluctuation effect even in Time Sequential Pro- 
cedure. These effects intermingle with each other 
and wc could not discuss them separately. However, 
it can be technically overcome in Time Sequential 
Procedure, while, in principle, it seems to be almost 
difficult to be solved in Vcut Procedure. 

Up to now, we restrict our discussion around 
electromagnetic cascade shower to Bethe-Heitler 
shower. However, we could not neglect LPM effects, 
related to the interpretation of extremely high en- 
ergy muon events in future. One is related to elec- 
trons [42-47] and other is related to muon [48,49]. 
We cannot neglect the LPM effect on the behav- 
iors of electromagnetic cascade themselves above 
- lO^^eV in Water [45]. However, the LPM effect 
is supposed to be effective above lO^^eV in the case 
of muon induced electromagnetic cascade showers. 

In such extremely high energies, range fluctu- 
ation of muon may alter their feature essentially 
compared with that of present situation. Further- 
more, above lO^^eV, we cannot neglect the LPM 
effect related to the muons [48,49]. Namely, above 
lO^^eV, Cherenkov light spectrum is supposed to 
become essentially different from those in the ab- 
sence of two kinds of LPM effects at present. 

In the present paper, we restrict our discussion to 
muons themselves in high energies. In subsequent 
papers, we will extend our discussion Cherenkov 
light yield via electromagnetic cascade showers from 
different interactions, such as, bremsstrahlung, di- 
rect electron pair production and photonuclear in- 
teraction, taking into account of muon energy spec- 
trum for imaging KM3 detector. Here, main sub- 
ject will be the discussion of Cherenkov light yeild 



from muon induced electromagnetic cascade show- 
ers. 

In discussion of the problems around high energy 
neutrino spectrum from the universe, it should be 
noticed that the reliable results are obtained only 
through the utilization of the stochastically correct 
tools, taking into account of very few number of 
experimental events, in addition to strong steepness 
of high energy neutrino spectrum. 
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